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Abstract 

We investigate non-standard Hamiltonian effects on neutrino oscillations, which are effective 
additional contributions to the vacuum or matter Hamiltonian. Since these effects can 
enter in either flavor or mass basis, we develop an understanding of the difference between 
these bases representing the underlying theoretical model. In particular, the simplest of 
these effects are classified as "pure" flavor or mass effects, where the appearance of such a 
"pure" effect can be quite plausible as a leading non-standard contribution from theoretical 
models. Compared to earlier studies investigating particular effects, we aim for a top- 
down classification of a possible "new physics" signature at future long-baseline neutrino 
oscillation precision experiments. We develop a general framework for such effects with 
two neutrino flavors and discuss the extension to three neutrino flavors, as well as we 
demonstrate the challenges for a neutrino factory to distinguish the theoretical origin of 
these effects with a numerical example. We find how the precision measurement of neutrino 
oscillation parameters can be altered by non-standard effects alone (not including non- 
standard interactions in the creation and detection processes) and that the non-standard 
effects on Hamiltonian level can be distinguished from other non-standard effects (such as 
neutrino decoherence and decay) if we consider specific imprint of the effects on the energy 
spectra of several different oscillation channels at a neutrino factory. 
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1 Introduction 



Neutrino physics has entered the era of precision measurements of the fundamental neutrino 
parameters such as neutrino mass squared differences and leptonic mixing parameters, and 
neutrino oscillations are the most credible candidate for describing neutrino flavor transi- 
tions. Nevertheless, there might be other sub-leading mechanisms participating in the total 
description of neutrino flavor transitions. Thus, in this paper, we will investigate such mech- 
anisms on a fundamental level, which will give rise to non-standard effects on the ordinary 
framework of neutrino oscillations. 

In a previous paper [1], we have studied non-standard effects on probability level based on 
"damping signatures" , which were phenomenologically introduced in the neutrino oscillation 
probabilities. However, in this paper, we will investigate so-called non-standard Hamiltonian 
effects, which are effects on Hamiltonian level rather than on probability level. Recently, 
three different main categories of non-standard Hamiltonian effects have been discussed in 
the literature. These categories are non-standard interactions (NSI), flavor changing neutral 
currents (FCNC), and mass varying neutrinos (MVN or MaVaNs). In addition, other effects 
which result in effective additions to the Hamiltonian have been studied, such as from extra 
dimensions [2]. Below, we will shortly review the categories of effects which can be studied 
using this framework. 

In general, in many models, neutrino masses come together with NSI, which means that the 
evolution of neutrinos passing through matter is modified by non-standard potentials due to 
coherent forward-scattering of NSI processes u a + / — > v@ + /, where a, (3 = e, //, r and / is 
a fermion in matter The effective NSI potentials are given by Vnsi = V^GpN^p, where 
Gp is the Fermi coupling constant, N d is the down quark number density, and e a ps are 
small parameters describing the NSI [3]. See, e.g., Ref. [4] for a recent review. Furthermore, 
matter-enhanced neutrino oscillations in presence of Z-induced FCNC have been studied in 
the literature [5-7]. See also, e.g., Refs. [8,9] for some earlier contributions. Especially, NSI 
and FCNC have been investigated in several references for many different scenarios such 
as for solar [10-13], atmospheric [14-18], supernova [19], and other astrophysical neutrinos 
as well as for CP violation [20], the LSND experiment [21], beam experiments [22], and 
neutrino factories [23-28]. 

The idea of MVN was proposed by Fardon et al. in Refs. [29,30]. This idea is based on the 
dark energy of the Universe being neutrinos which can act as a negative pressure fluid and 
be the origin of cosmic acceleration. Furthermore, several continuation works on MVN have 
been performed in the context of scenarios for the Sun and the solar neutrino deficit [31,32], 
but also in various other contexts [33-45] . In addition, it should be mentioned that neutrinos 
with variable masses have also been studied earlier than the idea of MVN [29,46-49]. 

While earlier studies have discussed individual theoretical models and their effects on future 
neutrino oscillation experiments (bottom- up), our approach will be the top-down. We start 
from general assumptions to investigate the properties of non-standard Hamiltonian effects, 

1 Note that, in general, the production and detection vertices could also be modified. However, in this 
paper, we focus on the neutrino oscillation probabilities which, in the limit of ultrarelativistic neutrinos, 
decouple from the creation and detection processes. 



1 



and later apply them to specific models and discuss how to identify individual effects. The 
goal of this approach is the classification of a possible "new physics" signature in future long- 
baseline neutrino oscillation experiments. Although it is very likely that such a signature 
will fit many different non-standard models, it has hardly been discussed in the literature 
how to distinguish (even qualitatively) different theoretical models which could all describe 
this effect, and what the methods for that identification could be. For this purpose, we make 
rather unspecific assumptions for the particular type of effect and rather assume that the 
theoretical model will predict a leading effect which can be considered to be of a "simple" 
form in a specific basis ("pure" effect), which can be either flavor (or mass) conserving or 
flavor (or mass) violating. 

The paper is organized as follows: First, in Sec. [2], we define non-standard Hamiltonian 
effects as effective additional contributions to the vacuum Hamiltonian similar to matter 
effects. The definition is performed for n neutrino flavors. Next, in Sec. [31 we specialize our 
discussion to two neutrino flavors, where we derive the effective neutrino parameters as well 
as the resonance conditions in both flavor and mass bases including non-standard Hamil- 
tonian effects. We also discuss experimental strategies to test and identify non-standard 
Hamiltonian effects at the example of v e «-> flavor transitions. Then, in Sec. HI we study 
some aspects of the generalization to three-flavor case, whereas in Sec. [5j we give a numer- 
ical example of how non-standard Hamiltonian effects can affect a realistic experimental 
setup and discuss how to tell non-standard Hamiltonian effects apart from damping effects. 
Finally, we summarize our results and present our conclusions in Sec. O 



2 Parameterization of non-standard Hamiltonian effects 



In the standard neutrino oscillation framework with n flavors, the Hamiltonian in vacuum 
is given by 

H = —Udmg(m 2 1 ,m 2 2 ,...,ml)U t (1) 

in flavor basis, where E is the neutrino energy, U is the leptonic mixing matrix, and m; 
is the mass of the ith neutrino mass eigenstate. Any Hermitian non-standard Hamiltonian 
effect will alter this vacuum Hamiltonian into an effective Hamiltonian 

H e s = H + H , (2) 

where H' is the effective addition to the vacuum Hamiltonian. We note that this reminds 
of neutrino mixing and oscillations in matter [8] with H' given by a diagonal matrix with 
the effective matter potentials on the diagonal, i.e., 

H' = H mat = di&g{V, 0, . . . , 0) - -±=G F N n l n , (3) 

where V = y/2GpN e is the ordinary matter potential, Gf is the Fermi coupling constant, 
N e is the electron number density (resulting from coherent forward-scattering of neutrinos), 
N n is the nucleon number density, and l n is the n x n unit matrix^. Just as the presence 

2 If sterile neutrinos are present, then there is no interaction between the sterile neutrinos and the matter 
through which they propagate. Thus, the 1„ term is replaced by a projection operator on the active neutrino 
states. 
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of matter affects the effective neutrino mixing parameters, the effective neutrino mixing 
parameters will be affected by any non-standard Hamiltonian effect. In the remainder of 
this text, we will treat the effective Hamiltonian 

H c s = H + H + H mat , (4) 

i.e., we will treat the non-standard effects along with the matter effects. However, in Sec. HJ 
we treat only the part Hq + H' in order to obtain the parameters of the Hamiltonian to which 
the standard matter effects are then added. Since standard matter effects are generally taken 
into account, H + H 1 will be mistaken for the vacuum Hamiltonian H if the non-standard 
effects are not considered. 

Since any part of the effective Hamiltonian that is proportional to the nxn unit matrix only 
contributes with an overall phase to the final neutrino state, it will not affect the neutrino 
oscillation probabilities. This means that we may assume H' to be traceless and also that 
we may subtract tr(H)/n from the effective Hamiltonian to make it traceless. Any traceless 
Hermitian nxn matrix A may be written as 

N 

A = ^QAi, (5) 

i=l 

where the q's are real numbers, the Aj's are the generators of the su(n) Lie algebra (i.e., 
A is an element of the Lie algebra), and N = n 2 — 1 is the number of generators. Hence, 
clearly, any non-standard Hamiltonian effect H' is parameterized by the n 2 — 1 numbers 
q. In summary, we choose the coefficients of the generators of the su(n) Lie algebra to 
parameterize any non-standard Hamiltonian effect. 

Furthermore, in any basis (e.g., flavor or mass basis), we may introduce su(ri) generators 
Aj such that n(n — I)/2 generators are off-diagonal with only two real non-zero entries, 
n(n — f)/2 generators are off-diagonal with only two imaginary non-zero entries, and n — 1 
generators are diagonal with real entries. For example, in the case of n = 2, we have the 
Pauli matrices 

Ax = or = (J I), A 2 = a 2 = ( ° J , A 3 = a 3 = ( J ^ J . (6) 

We will denote the set of generators which are of the form A, in flavor basis by pi and the 
set of generators which are of this form in mass basis by Tj. Obviously, in flavor basis, we 
have the relations 

Pi = Xi and n = U\iU ] , (7) 
where, in the case of two neutrino flavors, 

/ cos(#) sin(0) \ 
V - sin(0) cos(0) J 

is the two-flavor leptonic mixing matrix and 9 is the corresponding mixing angle (when 
treating the three-flavor case, we will use the standard parameterization of the leptonic 
mixing with three mixing angles $12, 023; 0i3 ; and one CP violating phase 5qp)- This implies 
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that pi and would be equal if there was no mixing in the leptonic sector. Furthermore, it 
is obvious that the matrices can be written as linear combinations of the matrices and 
vice versa. Therefore, there is, in principle, no difference between effects added in flavor or 
mass basis if one allows for the most general form of the non-standard contribution. 

We now define any non-standard effect as a "pure" flavor or mass effect if the corresponding 
effective contribution to the Hamiltonian is given by 



respectively, where c G M. This means that we restrict the "pure" effects to be of very 
specific types, where the actual forms are very simple in a given basisJl Given the possi- 
ble theoretical origin, this approach is quite plausible if one assumes that the underlying 
theoretical model will produce one leading flavor (or mass) changing or conserving effect. 
Generally, the parameter c can depend on many different quantities, e.g., the matter den- 
sity or the neutrino energy. In particular, the dependence on the neutrino energy ( "spectral 
dependence" ) may allow for the unambiguous identification, or, in the case of mass- varying 
neutrinos, the matter density dependence may indicate this type of effect. However, any 
approach investigating such dependencies has to use specific models, and the actual repre- 
sentation by Nature may easily be overseen. Therefore, we do not require this information 
in this study and rather investigate the generic impact of effects in the flavor or mass basis. 
In addition, we note that the matter density or energy dependence of the non-standard 
effects should be very weak for a given terrestrial experiment with a specific matter density 
profile. Only for effects motivated by MVN, i.e., mass effects, we will use the same energy 
dependence as for the masses themselves for numerical simulations. In general, if a large 
span of energies is available, one should of course also try to distinguish different specific 
models through their different energy dependencies. 

This choice of pure effects implies that only one of the generators of the Lie algebra is present, 
since a general linear combination, such as Eq. ([5]), can always be interpreted in both bases. 
Thus, we define a flavor or mass conserving (violating) effect as any effect where the effective 
contribution to the Hamiltonian is diagonal (off-diagonal) in the corresponding basisJl We 
note that a pure flavor (mass) violating effect corresponds to some interaction between 
two flavor (mass) eigenstates. For example, the su(2) generators p\ and pi correspond to 
flavor violating (or changing) effects, whereas p^ corresponds to flavor conserving effects. In 
summary, if we detect an arbitrary non-standard effect, it is the simple form in flavor or mass 
basis which makes it a flavor or mass effect by our definition. This approach can be justified 
by the fact that the simplest models for non-standard effects from the underlying theory 
correspond to specific patterns for the effective addition to the Hamiltonian. Therefore, our 
definition of a "pure" effect is a conceptually new one and it refers to a class of effects, which 

3 Because of our choice to use the Pauli matrices, a "pure" effect corresponds to the interaction of two 
flavor or mass eigenstates. This is also the reason for choosing to work with the Pauli matrices. In addition, 
it is also interesting to keep the real and complex parts of the off-diagonal entries separate (i.e., not working 
with the complex matrix elements directly, but rather a set of real parameters) in order to investigate the 
possibilities of probing CP violation effects. 

4 Strictly speaking, our definition distinguishes (in two-flavors) off-diagonal additions proportional to Xi 
(real) or A2 (complex). 
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can be interpreted in different ways. However, since simplicity is a basic concept in physics, 
this concept allows the choice of the most "natural" non-standard effects for further testing. 

The case of non-standard Hamiltonian effects on three-flavor neutrino oscillations, i.e., the 
case of n = 3, is quite similar to the one described above for the two-flavor case. Instead of 
the Pauli matrices, which are a basis of the su(2) Lie algebra, we now have to use the eight 
Gell-Mann matrices, which span the su(3) Lie algebra. Out of the Gell-Mann matrices, three 
are off-diagonal with two real entries, three are off-diagonal with two imaginary entries, and 
two are diagonal with real entries. Even though the principle of the three-flavor case is the 
same as that of the two-flavor case, it introduces many more parameters (more leptonic 
mixing angles, the complex phase in the leptonic mixing matrix, the extra mass squared 
difference, and the extra degrees of freedom for the non-standard effects), and therefore, 
turns out to be much more cumbersome to handle than the two-flavor case. In the following, 
we will therefore start by treating the two-flavor case in some detail and then continue by 
studying the similarities and differences when approaching the full three-flavor case. 

As far as the classification of current models in our notation is concerned, NSI and FCNC 
will be flavor effects, whereas MVN will produce mass effects. In general, NSI can be of two 
types: flavor changing (FC) and non- universal (NU) [4] . The off-diagonal elements of the 
effective NSI potential e al 3, where a ^ f3, correspond to FC, whereas the differences in the 
diagonal elements e aa correspond to NU. In addition, FCNC are flavor violating effects and 
MVN can be mass conserving. In principle, for our purposes, there is no difference between 
FC NSI and FCNC. 



3 Non-standard Hamiltonian effects in the two-flavor limit 

In this section, we study the general implications of non-standard Hamiltonian effects in 
the two-flavor limit. We discuss the effective parameter mapping including non-standard 
effects, and then, we apply it to a two-flavor limit as an example. 

3.1 Parameter mapping in two flavors 



In Appendix |A} we describe the general formalism of the two-flavor scenario, which can 
be used to obtain the results in this section. First, we discuss effects given in flavor basis, 
which are effects expanded in pi [cf., Eq. (jjTJ)]. In this case, flavor conserving effects will 
be contributions to the total Hamiltonian on the form H' = -F 3 p 3 , where F 3 G R, whereas 
flavor violating effects will be contributions on the form H' = Fipi + F 2 p2, where Fi G R. 
In flavor basis, the new effective parameters are given by 



Am 2 = Am 2 (, 

[® + M2*)] 8 + (£S0 ! 



sin' (26) 



where 



AEFa 

— 4+ sin 20) 



1 ( \F.F> N 2 



Am 2 



+ 



2VE AEF 3 

1 

Am 2 Am 2 



cos(2fl) 



(9) 
(10) 

(11) 
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is the normalized length of the Hamiltonian vector (see Appendix |A|) , Am 2 is the effective 
mass squared difference in flavor basis, and 9 is the effective mixing angle in flavor basisjf] 
In addition, the resonance condition is found to be 



2VE AEF 3 

1 

Am 2 Am 2 



cos(2fl), 



(12) 



which is clearly nothing but a somewhat modified version of the Mikheyev-Smirnov-Wolfen- 
stein (MSW) resonance condition [8,50,51]. From the resonance condition in Eq. (fT2l . it is 
easy to observe that the resonance is present for some energy E if and only if 



sgn(Am 2 ) sgn(V r ) sgn(l + 2F 3 /V) = sgn[cos(20)], 



(13) 



where sgn(Am 2 ) is dependent on the mass hierarchy, sgn(V^) is dependent on if we are 
studying neutrinos or anti-neutrinos, and sgn(l + 2F 3 /V) is dependent on the ratio between 
F 3 and the matter potential V [sgn(l + 2F 3 /V) being equal to —1 if and only if F 3 has 
a magnitude larger than \V/2\ and is of opposite sign to V}. Note that if there are flavor 
violating contributions added to the Hamiltonian, then these do not change the resonance 
condition. The sign of cos (28) can be made positive by reordering the mass eigenstates in 
the case of two neutrino flavors. However, we keep this term as it is, since this is not possible 
in the case of three neutrino flavors. This resonance condition can be easily understood, 
since the effective contribution to the Hamiltonian from any flavor violating effect will be 
parallel to the H 3 = plane, i.e., these contributions are off-diagonal. 

If we choose to describe the non-standard addition to the Hamiltonian in the mass eigenstate 
basis, then we find that the mixing parameters are given by 



Am 2 = Am 2 £, 

[ff cos(20) + (l-fM,)sin(20)] 2 +(f^) 2 



sin 2 (26) 



(14) 
(15) 



where 



2VE . „. AEM 1 
-sin(20) + 1 - 



Am 2 



Am 2 



1 ( 4EM 2 N 2 



+ 



Am 2 



+ 



2VE . „. AEM 3 



1 2 



and the resonance condition becomes 



2VE 4EM 1 



Arri 



+ 



Am 2 



sin(2#) + 



4EM 3 
Am 2 



cos(2#) = cos(2#). 



(16) 



(17) 



Note that both mass conserving effects and mass violating effects enter into the resonance 
condition, whereas only the flavor conserving effects entered in the corresponding expression 
in the flavor basis [c/., Eq. (jT2|) ]. This is due to the fact that the changes of the Hamiltonian 
vector from such effects are not parallel to the H 3 = plane (in flavor basis, see Appendix lAl). 
i.e., both of these effects affect the diagonal terms of the total Hamiltonian. However, M 2 
does not enter into the resonance condition, since a 2 = r 2 , i. e., the change of the Hamiltonian 
is off-diagonal also in the flavor basis. 



'Note that Fi may also change the effective Majorana phase. 
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3.2 Interpretation of experiments in the two-flavor limit 



Since a general analytic discussion of three-flavor neutrino oscillations including non-standard 
Hamiltonian effects would be very complicated, we focus on two neutrino flavors in this sec- 
tion. This approach can be justified if one assumes that the other contributions are exactly 
known or the two-flavor probabilities dominate. Of course, for short-term applications, 
small non-standard effects might be confused with other small effects such as sin 2 (26*13) 
effects [27]. Thus, a comprehensive quantitative discussion would be very complicated at 
present. 

In three-flavor neutrino oscillations, we can construct several interesting two-flavor limits 
of the probabilities P a p including non-standard effects related to two-flavor neutrino oscil- 
lations (see, e.g., Ref. [52]): 

ft. 1 --»'(*,) -n'P§£), (18, 

P. - o 1 -*,»(!*„) (19) 

P„ e — . s in 2 (29 13 ) sin 2 (^^) sin 2 (fe), (20) 

— l-sin 2 (2Msin 2 f^|^y (21) 

Note that all of these probabilities also contain the standard matter effects except from 
P^. In general, the su(3) generators (the Gell-Mann matrices) will give the degrees of 
freedom for non-standard Hamiltonian effects with three flavors. However, when studying 
the effective two-flavor neutrino oscillations, we only use the Gell-Mann matrices which are 
the equivalents of the Pauli matrices in the two-flavor sector that is studied. In addition, 
one can create two-flavor limits for oscillations into sterile neutrinos, such as in Ref. [2]. In 
the following, we will focus on small mixing and the case of Eq. fl20|) for illustration. We 
discuss the large mixing case in Appendix [Bl In addition, see Appendix [C] for subtleties 
with the definitions of the effective two-flavor scenarios. 

For small mixing, such as for Eq. (|20|) . we show in Fig. [1] the neutrino oscillation appearance 
probability P a $ for two flavors with small mixing, where the effects of the F^s are parame- 
terized relative to the matter effects (i.e., "1" on the vertical axis corresponds to an effect 
with Fi = V and "0" to no non-standard effects). In this figure, many of the following ana- 
lytic observations are visualized. The resonance condition in Eq. ffT2l) can always be fulfilled 
for the matter resonance (F3 = 0) by an appropriate choice of energy, baseline, neutrinos 
or antineutrinos, and oscillation channel. Obviously, we can read off from Eq. ffTUl) that 
at the resonance sin 2 (2(9) — > 1, where the matter resonance condition can be influenced by 
F 3 according to Eq. (|12|) . Therefore, the magnitude of sin 2 (2$) at the resonance (but not 
necessarily P a p) is independent of Fx, F2, and F3 by definition. However, F3 can shift the 
position of the resonance (such as in energy space). If we choose an energy far above the 
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Figure 1 : The two-flavor appearance probability P a p as a function of energy and the flavor 
conserving/violating fraction /j = Fj/V (normalized relative to matter effects). For the 
values of the neutrino parameters, we have used 9 = 0.16 ~ 9.2°, Am 2 = 0.0025 eV 2 , 
L = 3000 km, p = 3.5g/cm 3 , neutrinos only, Am 2 > 0, and > 0. 

resonance energy and Fj/V <C 1 (i = 1,2,3), then we have 



sin 2 (29) 



2VE _|_ 4EF 3 
Am 2 Am 2 



I]' 



(22) 



This means that Fx and F 2 can, for large enough energies, enhance a flavor transition, i.e., 
they increase the oscillation amplitude. In principle, one could distinguish F% from F 2 by 
a measurement at two different energies, because the mixed term from the square in the 
numerator of Eq. fl22l) has a linear (instead of quadratic) energy dependence. In practice, 
such a discrimination should be very hard. In addition, the quantity F3 can play the same 
role as the matter potential V, i.e., it can change the flavor transition for large energies. It 
is also obvious from Eqs. ffTUl) and (fTTl) that F 3 can affect the matter resonance energy and 
that it is directly correlated with the matter potential V, i.e., one cannot establish effects 
more precisely than the matter density uncertainty. 

In Sec. 13. 1^ we have also discussed mass effects, such as coming from MVN. Since a pure 
Mi or M 3 effect translates into a combination of F% and F 3 [cf., Eq. fl39l ], we expect to 
find a mixture of F\ and F3 effects, i.e., both F\ and F3 effects have to be present. Thus, 
if we assume that there is only one dominating "pure" non-standard contribution (F\, F 2 , 
F 3 , Mi, M 2 , or M 3 ), then this simultaneous presence points toward a mass effect. Clearly, 
an M 2 effect, on the other hand, cannot be distinguished from an F 2 effect [cf., Eq. (J39~l) ]. 
A different property of M 3 , which is not so obvious from Sec. 13.11 but very obvious already 
from Eqs. (pQ), ([2]), and ((Tj): Since M 3 is diagonal in mass basis, it corresponds to an energy 
dependent shift of the vacuum mass squared difference. As a consequence, in vacuum, the 
effective mixing angle is not modified by M 3 [cf., Eq. (1151) ]. Thus, the oscillation amplitudes 
are not modified by M 3 , but the oscillation pattern shifts (contrary to F 3 effects, where also 
the amplitude changes). In this case, the resonance condition becomes meaningless and 
the amplitude becomes sin 2 (29) = sin 2 (29). Note that a direct test using one experiment 
only makes it hard to identify mass effects uniquely if they are introduced with the same 
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energy dependence as the vacuum masses (because they can be rotated away by a different 
set of neutrino oscillation parameters). Thus, other methods might be preferable, such as 
modified MSW transitions in the Sun [31,32] or reactor experiments comparing air and 
matter oscillations [53]). 

Another class of effects has been discussed by Blennow et al. in Ref. [1]. In this study, so- 
called "damping effects" could describe modifications on probability level instead of Hamil- 
tonian level (such as neutrino decay, absorption, wave packet decoherence, oscillations into 
sterile neutrinos, quantum decoherence, averaging, etc.). It is obvious from Eq. (3) in Ref. [1] 
that these damping effects do not alter the oscillation frequency, while we can read off from 
Eqs. AH]) and (jllj) that it is a general feature of non-standard Hamiltonian effects that the 
oscillation frequency is changed. However, for damping effects, the oscillation amplitude 
can be damped either by a damping of the overall probability ( "decay-like damping" ) or by 
the oscillating terms only ( "decoherence-like damping"). In the first case, the total proba- 
bility of finding a neutrino in any neutrino state is damped for all energies, whereas in the 
second case, it is constantly equal to one while the individual neutrino oscillation probabil- 
ities are damped in the oscillation maxima and enhanced in the oscillation minima. Since 
all (small) effects one could imagine in quantum field theory, involving the modification of 
fundamental interactions or propagations, can be described by either coherent or incoherent 
addition of amplitudes, one can expect that the two classes of Hamiltonian and probability 
(damping) effects can cover all possible effects. However, in practice, potential energy, en- 
vironment, and explicit time dependencies (such as from a matter potential) can make life 
more complicated. 



4 Three-flavor effects 



As was stated in the Sec. [21 the general three-flavor case is quite complicated. However, if 
we assume that the non-standard effects are small, then we can use perturbation theory to 
derive expressions for the change in the neutrino oscillation parameters. For example, the 
elements of the effective mixing matrix are given by 

Uai = (v a \vi) , (23) 

where is the eigenstate of the full Hamiltonian. To first order in perturbation theory, 
we have 

M = l"*> + E %^ M - l"«> + 2E H^r h) , (24) 

... EJi J- J j ... t-*ii<>.ji 

and thus, we find 

H' 

U ai ~ U ai + 2E V (25) 
or, in terms of the non-standard addition given in flavor basis, 

U ai ^ ai + 2£VV U ^ ; f^ U aj . (26) 
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We note that this approach is valid only if \2EH-j/ Am 2 ^ <C 1. If this is not valid, then we 
have to use degenerate perturbation theory in order to obtain valid results. 

It was discussed in Refs. [26,27], that if is small enough, then possible NSI in the creation, 
propagation, and detection processes may mimic the effects of a larger 13 (this can also 
be the case for other effects which are not usually treated along with neutrino oscillations, 
such as damping effects [1]). Here, we consider only the propagation effects separately and 
consider how this alone could affect the determination of #13. The reason for doing so is that, 
while NSI can also affect the creation and detection processes, other non-standard effects, 
e.g., MVN, may not. With the perturbation theory approach described above, this becomes 
very transparent, and is probably one of the most interesting applications of non-standard 
effects. In any experimental setup, the value of the mixing angle #13 is determined by the 
modulus of the element U e3 of the neutrino mixing matrix U. If we include non-standard 
effects, then the effective counterpart of this element is given by 



IE 

U e3 ~ u e3 + -^ r (l + as 2 12 )(s 23 H' eijt + c 23 H' eT ) + 



31 



2E 



c 2 23 ^ T - s 2 23 H' Tfl + \ sin(2fl 23 )(^ - H*„ 



(27) 



where we have made a series expansion to first order in a = Am^/Am^ ~ 0.03 and 
disregarded terms of second order in both H' and #13. 

If U e3 is smaller than or of equal size to the other terms in this expression, then the #13 
determined by an experiment will not be the actual #13 unless the non-standard effects are 
taken into account. It is worth to notice that if 823 = 45°, then C23 = S23 and only the 
imaginary part of H' = (H')* will enter into the expression for U e3 , indicating that if 
the leading term is the one containing H' , then the effective CP violating phase will be 
±90°. Another interesting observation is that even if there are no non-standard effects, 
there is a term proportional to AV = H' — H' TT in this expression. Because of the different 
matter potentials for and v T due to loop-level effects, this quantity will be of the order 
AV ~ 10~ 5 V. 

In Fig. [21 we plot the possible range of \U e3 \ as a function of e max VE, where H' a/3 = e a pV , V 
is the matter potential, and \e a p\ < e max . For comparison, a neutrino factory with a neutrino 
energy of E = 50 GeV and a matter density of 3 g/cm 3 will have VE ~ 6 • 10~ 15 MeV 2 and 
the position at which we need to consider the possible range of \U e3 \ then depends on the 
bounds on the non-standard parameters e a p. In general, the bounds for e Q/ g depend on the 
type of non-standard effect and the types of interactions that are considered. In the case of 
NSI, it is common to write the non-standard interaction parameters as 



e 



a/3 



f 



where we sum over different types of fermions, e£g depends on the non-standard interaction 
with the fermion /, and Nf is the number density of the fermion /. In addition, e{g is often 
split into e£g = + e^, where L and R denotes the projector used in the fermion factor 
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le IVE [MeV 2 ] le IVE [MeV 2 ] 

max L J max L J 



Figure 2: The range of possible \U eS \ as a function of e iaax VE. The plots are arranged 
so that the left panels correspond to #13 = 10° and the right panels to #13 = 0, while the 
lower panels correspond to a non-standard effect with e er 7^ and the upper panels to a 
non-standard effect with e rr 7^ 0. The qualitative behavior for other values of #13 is similar 
to the behavior for 9 = 10°. (Note the different scales on the vertical axes.) 



of the effective non-standard Lagrangian density, i.e., 

C eS = -2V2G F J2 E €p(^lpL^)(h p Pf). (29) 

/ P=L,R 

Recent bounds for e^Jl can be found in Ref. [54] for electron neutrino interactions with 
electrons (i.e., e^) and Ref. [55] for interactions with first generation Standard Model 
fermions. As an example, the bounds from Ref. [54] for the e er (which is considered in 
Fig. ED are 

- 0.90 < ef T < 0.88 and -0.45 < < 0.44, (30) 

respectively. This means that the bounds, especially in this sector, are weak, which we will 
use in the next section. From Fig. [2j we can deduce that the off-diagonal e er terms have a 
larger potential of altering the value of \U e 3\ than the diagonal e TT terms, the maximal value 
can even exceed l/v2, corresponding to #13 = 45°. In addition, it is possible to suppress 
the effective #13 to zero if introducing non-standard effects. It follows that a relatively large 
9 1 3 signal, bounded only by the size of the non-standard effects, can be induced or that a 
large 6*13 signal can be suppressed by non-standard effects. Note that the effects quickly 
disappear at low energies, e.g., in reactor experiments. In order to tell a genuine # 13 signal 
apart from a signal induced by non-standard interactions, it is necessary to study the actual 
distortion of the energy spectrum induced by the neutrino oscillations. 
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5 A numerical example: Neutrino factory for large sin 2 (20i 3 ) 



This section is not supposed to be a complete study of non-standard Hamiltonian effects, 
but to demonstrate some of the qualitatively discussed properties from the last sections in a 
complete numerical simulation of a possible future experiment using the exact three-flavor 
probabilities. Therefore, we have to make a number of assumptions. We use a modified 
version of the GLoBES software [56] to include non-standard effects. As a future high- 
precision instrument, we choose the neutrino factory experiment setup from Refs. [57,58] 
with L = 3 000 km, a 50 kt magnetized iron calorimeter detector, 1.06 ■ 10 21 useful muon 
decays per year, and four years of running time in each polarity^ This experiment uses muon 
neutrino disappearance and electron to muon neutrino appearance as oscillation channels for 
both neutrinos and antineutrinos (in the muon and anti-muon operation modes combined). 
For the neutrino oscillation parameters, we use sin 2 2^12 = 0.83, sin 2 2^23 = 1, Am 2 x = 
8.2 • 10~ 5 eV 2 , and Am§ x = 2.5 • 10~ 3 eV 2 [59-62], as well as we assume a 5 % external 
measurement for Am^ and 9 12 [60] and include matter density uncertainties of the order 
of 5 % [63,64]. In order to test precision measurements of the non-standard effects, we use 
sin 2 (2# 13 ) = 0.1 close to the CHOOZ upper bouncf] [67], as well as we assume a normal 
mass hierarchy and Sqp = 0. For simplicity, we do not take the sgn(Am 2 1 )-degeneracy [68] 
into account, but we include the intrinsic (6*13, Sep) -degeneracy [69], whereas the octant 
degeneracy does not appear for maximal mixing [70]. Note that we do not include external 
bounds on the non-standard physics and sin 2 (2013), which, for instance, mean that we allow 
"fake" solutions of sin 2 (26 , 13 ) above the CHOOZ bound. This assumption is plausible, since, 
depending on the effect, the CHOOZ bound may have been affected by the non-standard 
effect as well. 



5.1 Test model 



Since we choose sin 2 (26*13) to be large, let us first of all focus on the appearance channel 
of u e oscillating into (or v e oscillating into v^). Expanding in small sin 2 (29 13) and a = 
Amli/Amjfl, we have for a — > (which should be a good approximation for sin 2 (26*13) ^ 
a 2 ~ 0.001) [71-73] 

o „ ,„ sin 2 [(l-i)A] . , 

P etl ~ sin 2 20 13 sin 2 #23—^ j-^, (31) 

(1 — A) 2 

where A = Am 2 , 1 L/(4E) and A = ±2\^2GFn e E/ Am^. Similarly, P er is described by this 
equation with sin 2 #23 replaced by cos 2 #23- This means that we may be effectively dealing 
with the two-flavor limits described in Sec. [3X depending on the degree the non-standard 



6 Compared to Ref. [57], we use a 2.5 % systematic normalization error for all channels as in Ref. [58]. 

7 In general, a large sin 2 (2#i3) will imply a large signal in the appearance channel. However, non-zero 
effective sin 2 (26'i3) could arise even if sin (2#i3) = 0, c/., Fig. [5] For effects which are diagonal in flavor 
basis, a large sin 2 (2#i3) would be preferred in order to make an observation of the non-standard effect. We 
have used large sin 2 (2# 13 ) as an example, since one may argue that the finding of new effects at present 
experiments (such as MINOS) may lead to a good reason for constructing a neutrino factory. One should 
also observe that, in principle, it would be possible to find non-standard effects at, e.g., MINOS [65,66]. 
However, the precision of a neutrino factory would be more sensitive to small effects, and thus, more useful 
for distinguishing between effects. 
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effects are different for the \i and r ffavors (c/., Appendix ICl). 

Using the parameterization in Eqs. ([5]) and ([6]) applied to the 1-3-sector, we therefore adopt 
the following Hamiltonian: 



H e f[ 




M 3 


Mi + LU- 



Ami 




21 



Mi - iM 2 


Am^ - M 3 



+ 



V + F 3 F 1 -iF 2 


F 1 +iF 2 -F 3 



(32) 



In this model, Mi and M 2 correspond to the CP conserving and CP violating parts of 
a mass-changing effect, whereas M 3 is a mass-conserving effect. In addition, F\ and F 2 
are the CP conserving and CP violating parts of a flavor-changing effect, whereas F 3 is a 
flavor- conserving effect. As motivated before, it is plausible to assume that one of these 
non-standard effects may be dominating the other ones, because many models predict such 
a dominating component and the experimental constraints on some quantities are rather 
strong. In addition, Eq. (1321) implies that the effects are mainly present in the 1-3-sector, 
which can be motivated by rather weak experimental bounds on the ^ r -sector. For example, 
the bounds on the matrix element H' eT are rather weak in the case of NSI, making it viable 
that this term is dominating the NSI Hamiltonian. In this case, we obtain 



Thus, we have a flavor violating effect with Fi representing the CP conserving part of the 
NSI and F 2 representing the CP violating part of the NSI. The form of the mass effects has 
been chosen to match the expected energy dependence of MVN in order to discuss effects 
with realistic spectral (energy) dependencies. 

Note that the parameterization in Eq. (1321) does not exactly correspond to the two-flavor 
limit even for a — > 0, since there are some non-trivial mixing effects in the 2-3-sector as 
described in Appendix O This parameterization is also obviously not the whole story in 
the three-flavor scenario. For instance, we assume the same sign for effects on neutrinos 
and antineutrinos, which may, depending on the model, not apply in general. However, we 
will demonstrate some of the characteristics from Sec. 13.21 with this approach. In addition, 
note that we have now adopted a specific energy dependence of the flavor and mass effects, 
where the definition of the energy dependence in the M's is slightly different from the one 
in the M's in Sec. [2J i.e., M = M/(2E). In this case, the mass effects could be coming 
from MVN changing the mass eigenstates, whereas the flavor effects correspond to some 
NSI approximately constant in the considered energy range. We will quantify the size of 
the Fi and Mj in terms of the normalized quantities fi = Fi/V (for p = 3.5g/cm 3 ) and 
Hi = Mi/Am 31 (for Am| x = 2.5 • 10 _3 eV 2 ). This quantification makes sense, since it is 
obvious from Eq. (I32p that the effect of these quantities will have to be compared with the 
order of V and Am^, respectively. Note that f\ — i/ 2 = £ e eT from Sec. HJ which means that 




V^Ime, 



(33) 
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Figure 3: Correlation between simulated models (columns) and fit models (rows). The 
areas of the disks represent the discovery potentials of the simulated "pure" effects (pa- 
rameterized in terms of fi or /ij) given that a different pure effect (fit model) is allowed 
(minimum value of a deviation from zero necessary in either direction for a 3a discovery). 
Therefore, the larger the disk, the more difficult it will be to distinguish a pure effect from 
another one. Note that we use cutoffs of < 0.3 and < 0.5 (largest gray disks), since 
some models cannot even be distinguished for much larger values. The areas of the rest of 
the disks are normalized with respect to these cutoffs for simulated flavor and mass effects. 



it will be interesting to compare the precisions of fi and / 2 to the current bounds for e er . 
Furthermore, note that the mass effects can be simply rotated away by a different choice of 
the mixing matrix and the mass squared differences because of the same energy dependence 
in this example. However, since we assume the solar parameters to be measured externally, 
we will observe that constraints to the Mj can be derived. Such an external measurement 
with a similar environment dependence to the neutrino factory comes from KamLAND, 
which turns out to be very consistent with the ones from solar neutrino experiments. Since 
most non-standard effects in oscillations are dependent on the matter density (such as 
MVNs with acceleron couplings to matter fields, or non-standard flavor- changing matter 
effects generated by higher-dimensional operators), it is plausible to assume that strong 
constraints hold for the solar sector because of the very different environments/densities 
within the Sun and the Earth. 

5.2 Identifying specific pure effects 

If we discover a non-standard effect, it will be an interesting question how easily it can 
be identified. Assuming one dominating effect of the mass or flavor type, which we have 
introduced as "pure effect", we want to know how well it can be distinguished from other 
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Quantity 


Lower limit (la) 


Upper limit (la) 


Lower limit (3a) 


Upper limit (3a) 


h 


-0.008 


0.008 


-0.025 


0.026 


h 


-0.003 


0.003 


-0.008 


0.008 


h 


-0.016 


0.016 


-0.049 


0.082 


/'i 


-0.176 


0.118 


-0.218 


0.211 




-0.105 


0.126 


-0.181 


0.212 


^3 


-0.015 


0.015 


-0.044 


0.090 



Table 1: Discovery limits for the parameters in Eq. ( 132]) as parameterized fi = Fi/V and 
/ij = Mj/Am^ from the neutrino factory simulation (including correlations). 

such effects of different qualitative nature. Therefore, in Fig. [31 we show the correlation 
between simulated and fit pure effects. For this figure, we simulate a pure effect (column) 
and fit it with a different one (row), i.e., we marginalize over the respective fi or /ij. The 
areas of the disks are proportional to the minimum simulated value necessary to establish 
a 3a effect, where we have chosen a cutoff of \fi\ < 0.3 and < 0.5 (corresponding to the 
largest gray disks) @ This means that the size of the disks measures the correlation between 
two pure effects and the ability to discriminate those. 

One can easily make a number of qualitative observations from Sec. 13 .21 quantitative. First, it 
is hard to discriminate between F\ and F 2 ( CP conserving and CP violating flavor-changing 
effects), since these effects are qualitatively similar and highly correlated with (9 13 (as we have 
tested). However, if Nature implemented a flavor- changing Fx or F 2 effect, then one could 
easily establish it against F 3 and the pure mass effects. In general, note that a discrimination 
between flavor and mass effects is rather easy because of their different spectral dependence 
in this example (such as between F 2 and M 2 ). The difference to F3 can be explained by 
the different flavor- conserving nature of F3. The results look somewhat different for the F3 
column: Because of the correlation with p and all of the neutrino oscillation parameters 
(see below), it will be hard to establish this effect. For the simulated mass effects, the scale 
is different, i.e., one cannot directly compare the M-columns with the F-columns. Again, 
the mass effects can be distinguished from the pure flavor effects to some extent. However, 
it is quite impossible to establish a mass effect against another one, since they can be easily 
simulated by a different set of mass squared differences and mixing parameters with the 
same energy dependence. The only reason why the pure mass effects can be established in 
this example at all is that we have imposed external constraints on the solar parameters as 
motivated above. 

5.3 Discovery of non-standard physics and potential for improvements 

A very important issue of any pure non-standard effect is its evidence compared to the stan- 
dard three-flavor oscillation scenario. Therefore, in Table [TJ we show the discovery reaches 
for the parameters from Eq. (132]) against the standard three-flavor neutrino oscillation sce- 
nario. This means that the shown pure effects are simulated and the standard three-flavor 

8 Note that, for instance, the gray disks for /1 and fi correspond to the order of magnitude of the upper 
bounds in Eq. ([50]) . which means that testing considerably larger effects does not make sense. 
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Fi F 2 F 3 




Figure 4: Main impact factors (impact greater than 5 %) for the test of specific simulated 
models (captions) against standard three-flavor neutrino oscillations (3a measurement). The 
neutrino oscillation parameters refer to correlations with the respective parameter, "Syst." 
refers to systematics, and "p" refers to the matter density uncertainty. The impact factors 
are defined as in Ref. [57] as relative improvement when the respective quantity is fixed 
(correlations) or systematics is switched off. 

neutrino oscillation parameters are marginalized over. Comparing the precisions of f\ and 
/2 with the numbers in Eq. (!30j) is impressive. However, these discovery reaches depend on 
sin 2 (26 , i 3 ) (and 5cp) and we have assumed a very large sin 2 (2#i3) = 0.1 (and #cp = 0). Note 
that the reach in f 2 is actually better than the one in /i, which is different from what is 
found in the two-flavor limit in Sec. 13. 21 The reasons are the mixing effects in the 2-3-sector 
and that F 2 is a non-trivial source of CP violation in the three-flavor case. 

Except from these sensitivities, which somewhat depend on the specific model, the behavior 
for neutrinos and antineutrinos, and so on, it may be of some interest to obtain hints how 
these reaches can be improved. In order to study this aspect, we show the so-called "impact 
factors" for the test of specific simulated models against standard three-flavor neutrino 
oscillations in Fig. |H These impact factors test the relative impact of the measurement errors 
on the neutrino oscillation parameters and systematics. In order to compute them, the non- 
standard discovery limits are evaluated with all neutrino oscillation parameters marginalized 
over, matter density uncertainties included, and systematics switched on (standard). In 
addition, in order to test a specific impact factor, one neutrino oscillation parameter is 
fixed at one time (or systematics is switched off), and the corresponding discovery reach for 
the non-standard effect is compared to the discovery reach including all uncertainties and 
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systematics. The difference between these to discovery reaches describes the impact of a 
particular measurement error (or systematics), and the relative impact in Fig. H] quantifies 
what one needs to optimize for in order to improve the discovery reach. For example, for 
M 3 (lower right pie), the error on Am^ is the main impact factor in our model, which needs 
to be improved to increase the M3 discovery reach. 

Again, a number of aspects from Sec. 13.21 can be verified. For F\ and F 2 effects, systematics 
is the main impact factor, since these flavor effects determine the overall height of the 
appearance signal and are not introduced with a specific spectral dependence (remember 
that we use a conservative overall normalization error of 2.5 %). For F 3 effects, we have 
earlier determined the matter density uncertainty as an important constraint. However, 
improving the knowledge on Am^, #13, or <5cp does have a similar effect, since the extraction 
of the individual parameters becomes easier. For the mass effects, we encounter a completely 
different behavior. Remember that we have defined the mass effects with the same energy 
dependence as the mass squared differences, which means that particularly M 3 is easily 
mixed up with Am^. On the other hand, Mi and M 2 are related to a flavor change in the 
appearance channel via the mixing matrix, i.e., the leptonic mixing angle 813. Therefore, 
it is not surprising that such a flavor change can be interpreted as either a mixing or a 
mass-changing effect. Compared to Sec. 13.21 there are also a number of differences coming 
from the three-flavor treatment (solar and CP effects) and the mixing in the 2-3-sector. 
These effects introduce additional correlations with #23 and Sep- However, they are also 
the reason why, for example, M 3 can be constrained at all from this experiment alone [in 
the pure two-flavor case or without external constraints on the solar parameters, it would 
be impossible to distinguish between a non-vanishing M3 and a different Am^ if the mass 
effects had the energy dependence assumed in Eq. fl32l ]. 

5.4 Comparison to damping effects 

In the context of the non-standard effect identification, a more general question is the 
ability to distinguish Hamiltonian effects and effects on probability level. The probability 
level effects lead to damping of the neutrino oscillation probabilities ( "damping effects" ) and 
were studied in detail in Ref. [1]. They may originate from decoherence, neutrino decay, or 
other physics mechanisms. In this section, we address this identification in somewhat more 
detail in a qualitative manner. A relatively new ingredient for this identification is the use of 
the "Silver" [y e — > u T ) channel at a neutrino factory [74,75]. It has been noticed [28,65,66] 
that the Silver channel probability can be greatly enhanced for non-standard Hamiltonian 
effects. This corresponds to what we have found in Sec. 13.21 i.e., the Silver channel, which is 
similar to the "Golden" [v e —>■ v^) channel when there are no non-standard effects, behaves 
as our two-flavor limit in Sec. 13.21 for large energies. 

In Fig. [5j we show the impact of different types of effects on the neutrino oscillation proba- 
bilities in the Golden channel P efJi , the disappearance channel -P^JE and the Silver channel 
P eT (shown in columns) at a possible future neutrino factory (relevant energy range shown). 
The different rows correspond to scenarios with F\ (flavor-changing without CP violation, 

9 The probability is actually the survival probability. However, this is the relevant probability 
when searching for disappearance rather than the disappearance probability 1 — P MM . 
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Figure 5: The impact of four different non-standard effects (in rows) on three different 
oscillation channels (in columns): Golden u e — * Disappearance — > and Silver 
v e — > v r . The different rows correspond to F\ (flavor- changing without CP violation, 
Hamiltonian-level) , F3 (flavor-conserving, Hamiltonian- level), decoherence, and neutrino 
decay. The different model parameters for the different curves are given in the individual 
plots, where the thick curves correspond to the standard neutrino oscillation scenario. 
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Hamiltonian-level) , F3 (flavor-conserving, Hamiltonian-level) , decoherence, and neutrino de- 
cay, respectively. The different model parameters for the different curves are given in the 
individual plots, where the thick curves correspond to the standard neutrino oscillation 
scenario. For a description of the decoherence and decay models, see Ref. [1]. In short, 
the decoherence model corresponds to standard wave packet decoherence, whereas the de- 
cay model assumes equal decay rates for all mass eigenstates (which can, for instance, be 
motivated by a degenerate mass spectrum). Note that this figure is shown for neutrinos 
only and that the comparison between the neutrino and antineutrino behavior can provide 
information on the underlying physics as well. However, this behavior is model-dependent. 

Figure is very useful to study the characteristics of different effects and to illustrate how 
the information from different neutrino oscillation channels can be used to disentangle them. 
First, it is important to note that it is difficult to construct a damping effect without large 
impact on the disappearance channel. Since the event rates in this channel are very high, 
it is probably the first place to look for non-standard physics. In addition, damping effects 
tend to suppress the Golden and Silver channel probabilities around the oscillation peak, 
which can, depending on the model, be very different for Hamiltonian-level effects. However, 
as it can be read off from Fig. [5j Hamiltonian-level effects produce, similar to sin 2 (2^13), 
the largest effect in the appearance channels!*"! In particular, the flavor- conserving effect 
Fi may enhance the silver probability as demonstrated in the two-flavor limit in Sec. 13.21 
Comparing all panels in Fig. 0, we expect that the combination of all channels serves as a 
good model discriminator because each of the shown models has a unique signature if all 
these channels are combined. For example, we have tested that adding a 5 kt OPERA-like 
emulsion cloud chamber for the silver channel at the same baseline as the golden channel 
improves the Fi discovery reach by about 60 % (for simulation details, see Ref. [76]) because 
of the silver channel signature at large energies. Therefore, we believe that this combination 
of different channels in combination with precise oscillation parameter measurements and 
spectral signatures can reveal non-standard physics. 

6 Summary and conclusions 

For future long-baseline neutrino oscillation precision measurements, such as neutrino fac- 
tories, it will be an important question how to identify a non-standard effect. While it is 
very likely that many theoretical models will fit such a "new physics" discovery, the clas- 
sification of models corresponding to this discovery from a phenomenological point of view 
will be very important for the planning of the following generation of experiments. Hence, 
there is strong interest in a top-down approach to non-standard physics tests, since the 
impact of future measurements on theory has to be assessed to promote the experiment. So 
far, mainly the bottom-up approach has been used, which is testing specific models in an 
experiment. Therefore, it has been one of the main goals of this work to demonstrate the 
identification and separation of individual phenomenological classes by generic arguments. 

10 The weak influence of the non-standard Hamiltonian level effects on the disappearance probability 
1 — Pnn is purely due to the fact that Eq. (I32p has been assumed for the non-standard Hamiltonian, where 
we have assumed the states to be unaffected. However, as mentioned earlier, there are also stronger 
bounds on any NSI involving z/ M . 
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In summary, we have studied non-standard effects on neutrino oscillations on the Hamilto- 
nian level. We have parameterized these effects in terms of the generators of the Lie algebra, 
and we have introduced them in flavor (such as coming from FCNC or NSI) and mass (such 
as coming from MVN) bases. As a trivial fact, there is, in principle, no mathematical dif- 
ference between these effects if one allows the most general form in each basis. Given the 
detection of a general non-standard effect on Hamiltonian level, it is therefore not possible 
to classify it as a flavor or mass effect without further assumptions or knowledge and, from 
an empirical point of view, the classification is a matter of definition. Therefore, we have 
defined "pure" effects as effects which are proportional to specific individual generators. 
Those correspond to pure flavor/mass conserving/violating effects, i.e., effects which affect 
particular flavor or mass eigenstates. This definition makes sense if one assumes that the 
underlying theoretical model causes one dominating non-standard effect. It is then the sim- 
plicity of the form in the respective basis which defines the effect to be of flavor or mass 
type. Therefore, the concept of these pure effects allows the choice of the most "natural" 
class of models for further testing, which is most appealing from the physics point of view. 

From the analytical point of view, we have studied the effects in the two-flavor limit. We 
have derived the modified mass squared differences and mixing angles (parameter mappings) 
as well as the modified resonance conditions including standard matter effects. In addition, 
we have discussed the application of this two-flavor limit to experiments, in particular, to 
the neutrino oscillation probability P eM . This probability can be described to a first approx- 
imation by a two-flavor limit for large sin 2 (2# 13 ), where the sin 2 (26 l 13 )-term dominates the 
CP effects. In addition, non-standard effects in the 1-3-sector have so far very poor limits 
(such as e er ) and the driving parameter sin 2 (2#i3) is unknown, which means that there is 
room for confusion between 9 13 and non-standard effects (see, e.g., Ref. [26,27]). We have 
found that there are several generic features for different types of effects. While any flavor 
violating pure effect can obviously change the transition probabilities, it does not affect the 
resonance condition/energy. However, a flavor conserving pure effect changes the resonance 
condition similar to matter effects and is highly correlated with the matter density. In addi- 
tion, it can suppress the flavor transition for large energies similar to matter effects - even in 
vacuum. Pure mass effects behave, in principle, as rotations of the flavor effects by the mix- 
ing angles, i.e., a pure mass effect will be observed as a linear combination of flavor effects. 
However, for a pure mass conserving effect, these flavor effects combine with special char- 
acteristics, since the mass effect is similar to an (energy dependent) change of the vacuum 
mass squared difference, i.e., it basically squeezes or stretches the oscillation pattern. Since 
in quantum field theory any non-standard effect may originate in the coherent (Hamiltonian 
effect) or incoherent ("damping" effect) summation of amplitudes, we have compared the 
non-standard Hamiltonian effects to the previously studied "damping" effects on probability 
level. We have found that these two classes can be distinguished by typical characteristics. 
Non-standard Hamiltonian effects shift the oscillation pattern, while "damping" effects, in 
general, do not. In principle, the different classes of non-standard Hamiltonian effects can 
be identified by their modification of oscillation amplitudes for large energies, the shift of 
the matter resonance, the comparison of different L/ E-ranges, etc. 

We have also studied some aspects of the three-flavor generalization of general non-standard 
Hamiltonian effects using perturbation theory as well as numeric calculations. By assuming 
small non-standard Hamiltonian effects, we have derived expressions for the effective matrix 
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elements using perturbation theory and observed how the confusion theorem between 9%s 
and non-standard effects described in Ref. [27] arises at the Hamiltonian level. Our numeric 
calculations show that non-standard effects can alter the determination of #13 significantly 
at higher energies, while still preserving a high accuracy at lower energies (cf., Fig. [2]). 

Eventually, we have demonstrated, at a numerical example for a neutrino factory, that many 
of these features can be found in a realistic experimental simulation using three flavors and 
specific spectral (energy) dependencies for the non-standard effects. For example, while it 
is simple to distinguish a flavor changing effect from flavor conserving or mass effects in 
general, mass effects are hard to establish as long as the neutrino oscillation parameters are 
not known from an independent source (such as with a different matter density for MVN). 
In addition, we have compared the obtainable discovery reaches for e er to the current limits, 
and we have found at least an order of magnitude improvement for large sin 2 (2^13) and 
5qp = 0. We have also compared the Hamiltonian level effects to damping effects and found 
that they can be distinguished by their specific alteration of the spectra in different neutrino 
oscillation channels. 

Since the past has told us that neutrinos are good for surprises, the high precision measure- 
ments at future neutrino oscillation experiments might as well reveal a detection of "new 
physics" beyond the Standard Model (extended to include massive neutrinos). Therefore, 
we conclude that one should include general strategies to look for non-standard effects in fu- 
ture neutrino oscillation experiments, where we have followed a top-down approach: Instead 
of testing particular models (bottom-up), we have assumed that some inconsistency will be 
found first. Secondly, one may want to classify this inconsistency to be either a Hamiltonian 
or a probability level ("damping") effect. Finally, individual models are identified which fit 
this effect. Since we do not know exactly what we are looking for, such an approach might 
be a clever search strategy, and it can be useful to promote an experiment as a discriminator 
among different classes of theoretical models. Future studies should demonstrate how such 
an approach can be most efficiently extended to three neutrino flavors, which neutrino oscil- 
lation channels are most suitable, and what the correlations with the existing fundamental 
neutrino oscillation parameters imply. 
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A General formalism for the two-flavor scenario 



Any two-flavor Hamiltonian can be written in flavor basis on the form 

H = H a, (34) 

where H 6 M 3 and a = (a±, o 2 , 03) is the vector of the three Pauli matrices (cf., the 
pictorial description of two- flavor neutrino oscillations in Ref. [77]). For a time- independent 
Hamiltonian, the time evolution operator is given by 

S(t) = exp(-iHt). (35) 

Using the relation (A ■ a) 2 = \ A\ 2 , one obtains 

S(t) = 1 2 cos(\H\t) - x^- sin(|Jf \t), (36) 

H J 

where 1 2 is the 2x2 unit matrix. This gives the two-flavor neutrino oscillation probabilities 
of the form 

p aa = tr[P + S(t)\ = 1 - sin 2 (20) sin 2 (H), (37) 
P a(3 = tr[PS(t)] = sin 2 (2fl) sin 2 (H), (38) 



where 



sin 2 (20) = 4^, ^d k = \H\ 
2 Lrz 



\ 



Here the H^s are the components of the Hamiltonian and 9 is the effective mixing angle. 

In the standard two-flavor neutrino oscillation scenario, Hi = sin(26 l )Am 2 /(4i?), H 2 = 0, 
and i?3 = V/2 — cos(26) Am 2 / (AE) . In general, the resonance condition, i.e., the condition 
for maximal effective mixing, is H% = 0. The Hamiltonian is represented as a vector in 
M. 3 , the third direction being the "flavor" eigendirection. The mixing is given by the angle 
between the Hamiltonian vector and the flavor eigendirection. The mixing is maximal, i.e., 
sin 2 (29) = 1, when the Hamiltonian vector is orthogonal to the flavor eigendirection, which, 
as expected, is equivalent to the resonance condition. 

The flavor and mass bases, and thus, the flavor and mass effects, are intimately associ- 
ated with each other. For the case of n = 2, i.e., for two neutrino flavors, any effective 
contribution to the Hamiltonian can be written in either flavor or mass basis, i.e., as 
H' = F x p\ + F 2 p 2 + F 3 p 3 or H' = Min + M 2 r 2 + M 3 r 3 . Since the effect must be the 
same regardless of the basis it is expressed in, we obtain the relations 

Fx = Mi cos(2#) - M 3 sin(2#) 

F 2 = M 2 (39) 



F 2 = M 2 

F 3 = Mi sin(2#) + M 3 cos(2#) 



from Eq. 0, i.e., one obtains Fi and _F 3 by rotating Mi and M 3 by the angle — 29 as well 
as one has F 2 = M 2 . Thus, the transformation in Eq. fl39l) relates flavor and mass effects 
and shows that they are linear combinations of each other. 
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B Non-standard effects for large mixing 



In this appendix, we concentrate on pure effects in the limit of large mixing. When the 
mixing goes to maximal, we have cos(2#) — > 0. This means that the resonance condition 
in Eq. (fT2"l) cannot be fulfilled for F 3 = (i.e., "matter resonance") at reasonably large 
energies Q From Eq. (fTUj) . we can easily observe that Fx and F 2 will not modify sin 2 (2$) at 
all in the absence of matter (and _F 3 ) effects (for example, in P w or a vacuum probability). 
Independent of matter effects, F 3 can increase the suppression of sin 2 (26) for large energies. 
If the resonance condition in Eq. (fl2l) is fulfilled, then sin 2 (26 l ) will be independent of Fx 
and F 2 . However, in the presence of matter effects (such as for P ee in the limit #13 — > 0), Fx 
and F 2 can reduce the matter effect suppression of sin 2 (26* ) for large energies (cf., Fig. [6]), 
i.e., they can increase the effective mixing. Eventually, it is obvious from Eqs. (Q and ffTTl) 
that the oscillation frequency is always increased for positive Fj/Am 2 's. 

It can also be interesting (and quite illuminating) to study how different pure (flavor or mass) 
effects affect the effective neutrino mixing and oscillations. In Fig. [6j we plot the effective 
mixing resulting from "pure" flavor and mass effects. From this figure, some features become 
quite apparent, e.g., some generic features are the shift in the resonance energy for F3, M\, 
and M3, the non-zero high-energy mixing for all effects but the flavor conserving effect F%, 
and the appearance of an anti- resonance - where sin 2 (2$) goes to zero for some finite energy 
- for F u Aft, and M 3 . 

The shift in the resonance energy is simply due to the shift in the if3-component of the total 
effective Hamiltonian in flavor basis (as was mentioned earlier, the resonance condition is 
H 3 = 0). The fact that there is no shift of the resonance condition for Fx and F 2 was also 
discussed earlier. 

The reason why the effective high-energy mixing generally turns out to be non-zero is also 
quite easy to realize. At high energies, the effective matter potential, which is diagonal in 
flavor basis, dominates over the vacuum Hamiltonian. As a result, the effective mixing is 
usually zero at high energies. However, if there is a non-standard effect with a correspond- 
ing effective addition to the Hamiltonian which is non-diagonal and is either constant or 
increasing with energy, then the effective mixing at high energies will be fully determined 
by the ratio of the non-standard effect and the matter potential. 

The anti-resonance appears when Hi = H 2 = in flavor basis. Since H 2 = in the standard 
neutrino oscillation scenario, it is apparent that this anti-resonance will occur for some value 
of Fx. In addition, since Afi and M3 are linear combinations of Fx and F3, the anti- resonance 
will also appear for Aft and Af 3 effects, as can be seen from the plots in Fig. O 

There are also some interesting features that are specific for different effects. First, for Fx 
effects, we note that the resonance condition is unchanged and that the mixing is constant as 
a function of energy for Fx/V = — tan(2#) /2 (the reason for this is that the sum of the non- 
standard Hamiltonian and the matter potential is proportional to the vacuum Hamiltonian). 
Then, for flavor conserving F3 effects, we note that these correspond to changes in the 

n Note that, in this appendix, we assume that matter effects determine the resonance energy and the 
non-standard effects are sub-leading contributions, which may shift the resonance energy. Thus, we refer to 
the "matter resonance" as the resonance condition in Eq. (|12[) for F3 = 0. 
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Figure 6: The effective mixing sin 2 (2$) as a function of 2VE/Am 2 and the ratio between 
the pure flavor (Fx and F3) or mass (Mi and M3) effects and the matter potential V. The 
horizontal lines correspond to no non-standard effect and the mixing along them is therefore 
the same in all panels. The vacuum mixing is assumed to be 6 = 30°, which is close to the 
present value of the solar mixing angle # 12 (see, e.g., Ref. [62]). See the main text for a 
more detailed discussion. 



effective matter potential. For F$/V < —1/2, we obtain an effective matter potential which 
is negative, resulting in a disappearance of the resonance. Next, for Mi effects, the mixing 
is constant when M\jV = — sin(2#)/2, in analogy with the Fx effects (again, the reason is 
that the sum of the non-standard Hamiltonian and the matter potential is proportional to 
the vacuum Hamiltonian). Also in analogy with the Fx effects is that there is a value of 
2VE/Am 2 , where the mixing does not depend on the Mx/V ratio. However, in the case 
of Mx effects, this is not the resonance mixing, but rather a mixing of sin 2 (26) = cos 2 (26), 
which appears at 2VE/Am 2 = 1/ cos(26 l ). Finally, for mass conserving M 3 effects, we note 
that the resonance disappears when 2cos(26)M^/V < —1. 

The reason why the equivalent F 2 and M 2 effects are not included is that these effects always 
lead to an increase in the effective mixing angle for all energies, and thus, those plots do 
not show as many interesting features as the plots included. In addition, we note that if the 
non-standard effects are energy dependent, then the effective mixing will be given by the 
mixing along some non-constant function of 2VE/ Am 2 in Fig. [6l 
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C Two-flavor limits of three flavor scenarios 



In this appendix, we discuss subtleties with the definition of the effective two-flavor scenarios 
introduced in Sec. 13.21 Remember that the effective two-flavor neutrino oscillation scenarios 
should be defined in terms of the effective two-flavor sector in question. For example, 
in the limit when Am^ — > 0, the effective two-flavor sector is spanned by v e and u a = 
~S23^n + c 23 v t- Thus, the limit can be considered as an exact pure two-flavor scenario only 
if the non-standard effects preserve the two-flavor limit (i.e., no off-diagonal terms mixing v e 
and v a with the remaining neutrino state vj, = c^v^ + S2z v t)- If the non-standard addition 
to the Hamiltonian is given by 

Kp = z*pV, (40) 
then the corresponding addition in the basis spanned by \y e , u b , u a } is 




H' = V\ c 23 e* T - s 23 el T A B J . (41) 



where 



A = c 23 e W* + S 23 e TT - S 23 C23(£ht + £*) 



7'T 

'23 c Atr ~r ^23<-23V c W ~~ € rr) , 
C — s 23 e uu + c\n€ TT + S23C23(e„ r + 6* 



B = c 23 e Mr - S2 3 e* + s 2 3C2 3 (e MM - e 7 



, 23 C MM ^ ^23^-tt * ^23>-23V c /iT "r c^;. 

From this relation, we deduce that the limit will be a pure two-flavor case if e a/ g = for 
all non-standard effects which do not involve v e and C23e eM — S23e eT = (which could be 
implemented by, e.g., 6*23 = 45° and e eA1 = e er ). In general, some of the conclusions for the 
two-flavor case will therefore not apply to three flavors. We have, in the numerical example 
in Sec. |5l demonstrated which of the conclusions that do hold. The case when 6 13 — ► is 
similar to the case described above, with the exception that the effective two-flavor sector 
is now spanned by u e and i/& instead of v e and v a . For the limit #13 — > and Am^ — > 0, 
there is no subtlety and the two-flavor sector is spanned by and v T . 
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